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Abstract 

We  consider  Maxwell  systems  with  instantaneous  conductivity  and  nonlocal  in  time 
(hysteretic)  polarization  laws  which  are  characteristic  of  dispersive  dielectric  media. 
We  formulate  such  systems  in  an  operator  theoretic  framework  and  show  that,  under 
certain  conditions  on  the  dielectric  response  function,  the  resulting  systems  generate 
Co  semigroups.  It  is  shown  that  multiple  Debye  polarization  models  are  included  in 
those  for  which  a  semigroup  formulation  is  possible. 


1  Introduction 


In  a  forthcoming  monograph  [BBLOO]  we  have  developed  a  theoretical  and  computa¬ 
tional  framework  for  electromagnetic  interrogation  of  dispersive  dielectric  media.  In 
that  work  we  show  that  one  can  take  a  time  domain  variational  or  weak  formulation 
of  Maxwell’s  equations  in  dispersive  materials  and,  in  the  context  of  inverse  prob¬ 
lems,  use  partially  reflected  polarized  microwave  pulses  to  determine  both  dielectric 
material  properties  and  geometry  of  bodies  (specifically  for  plane  waves  inpinging 
on  slab  geometries  in  paradyms  which  approximate  far  field  interrogation).  This  is 
done  in  configurations  involving  either  supraconductive  reflecting  back  boundaries  or 
acoustically  generated  virtual  reflectors. 

The  propagation  and  reflection  of  electromagnetic  waves  in  dispersive  dielectric 
media  is,  of  itself,  an  interesting  topic  of  investigation.  As  we  point  out  in  the  next 
section  (and  demonstrate  computationally  in  [BBLOO]),  the  underlying  dynamical 
systems  are  not  typical  of  either  standard  parabolic  or  standard  hyperbolic  (even 
with  the  usual  dissipation)  systems  and  are  hence  of  mathematical  interest.  In  this 
short  note,  we  consider  the  Maxwell  system  for  rather  general  dispersive  dielectric 
media  and  show  that  such  systems,  under  appropriate  conditions  on  the  polarization 
law,  generate  Co  semigroup  solutions.  These  results  are  presented  in  the  context 
of  the  1-dimensional  interrogating  systems  developed  in  detail  in  [BBLOO]  and  we 
invite  interested  readers  to  consult  that  reference  for  more  detailed  discussions  and 
development  of  the  underlying  model  employed  here. 
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2  Modeling  of  Dispersiveness  in  Dielectric  Media 


We  begin  with  time  domain  Maxwell’s  equations  in  second  order  form  (e.g.,  see 
[BBLOO])  for  the  electric  field  E  =  E(t.  z )  of  1-dimensional  polarized  waves 

E  +  —P+  -jc-  c2E"  =  ~—js  (2.1) 

Co  Co  Co 

where  c  =  -j==  is  the  speed  of  light  in  vacuum,  Jc  is  the  conduction  current  density, 
Js  is  a  source  current  density  and  P  is  the  electric  polarization  of  the  dielectric 
medium.  We  assume  very  general  constitutive  material  laws  for  the  polarization  and 
conductivity  given  by 

P(t,  z)  =  (gp  *  E)(t,  z)  =  J  gp(t  -  s,  z)E(s ,  z)dz  (2.2) 

Jc(t,z)  =  (gc*  E)(t,  z)  =  [  gc(t  —  s,  z)E(s,  z)dz  (2.3) 

Jo 

where  we  have  tacitly  assumed  that  E(t ,  z)  =  0  for  t  <  0  and  that  both  gp(£,  z)  and 
gc(^,z)  vanish  for  £  <  0.  With  these  assumptions,  the  integrals  in  equations  (2.2), 
(2.3)  are  equivalent  to  integration  over  all  of  (— oc,  oc)  and  thus  are  indeed  convo¬ 
lutions.  The  displacement  susceptibility  kernel  gp  (also  referred  to  as  the  dielectric 
response  function (DRF))  and  the  conductivity  susceptibility  kernel  gc  introduce  non¬ 
locality  in  time  in  the  polarization  and  conductivity  relationships  [APM89],  [Jac75] 
which  is  equivalent  to  frequency  dependence  of  the  dielctric  permittivity  e  and  con¬ 
ductivity  a  when  using  a  frequency  domain  approach.  We  assume  that  either  P  or 
Jc  or  both  may  contain  instantaneous  (local  in  time)  components  by  introduction  of 
5  distributions  in  the  kernels  gp  and/or  gc,  respectively. 

A  medium  is  dispersive  if  the  phase  velocity  of  plane  waves  propagating  through  it 
depends  on  the  frequency  of  the  waves  [Jac75,  Chap. 7],  [Ch89,  Chap. 8].  To  determine 
the  dispersive  nature  of  a  medium  described  by  equations  (2.1)-(2.3)  we  seek  plane- 
wave  solutions  of  the  homogeneous  analogue  of  (2.1)  of  the  form  E(t,  z)  =  EQe~'l(ut±KZ'> 
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which  travel  in  the  z  direction  and  have  wavelength  A  =  2'k/k.  The  phase  velocity 
vp  of  these  waves  is  the  speed  at  which  planes  of  constant  phase  move  through  the 
medium.  In  this  case  the  argument  cot  —  n z  is  constant  and 

dz  , 

Vp=dt=  (2-4) 

Seeking  plane  wave  solutions  of  the  form  in  (2.1)  is  equivalent  to  seeking 

solutions  of  the  form  E0e±ZKZ  in  the  frequency  domain  version  of  (2.1).  Thus  we  use 
the  Fourier  transform  in  (2.1)  and  obtain 

-u?E  -  —P  -  —Jc  -  c2E"  =  0  (2.5) 

Co  Co 

where  we  have  ignored  the  source  term  Js  and  where  the  overhat  will  represent  the 
Fourier  transform  throughout.  Since  we  see  from  (2.2)  and  (2.3)  that  P  =  gpE  and 
Jc  =  gcE ,  this  can  be  written 

c2E"  +  u2(1  +  ^  +  ^)E  =  0.  (2.6) 

U6q  Co 

We  note  that  (2.6)  is  the  generalized  Helmholtz  equation  [Jac75,  p.  271] 

E"  +  k2E  =  0  (2.7) 


with 


K 


2 


<1  +  AffcH  + 

c2  V  ue0  e0 


(2.8) 


which  has  solutions  E(u,z)  =  E0e±ZK^z.  It  follows  that  the  corresponding  time 
domain  solutions  are  our  desired  solutions  of  the  form  E(t ,  z)  =  Eoe~l^t±ZKZ^  where 
the  wavenumber  n  =  n(u)  will  in  general  depend  on  the  frequency  uj.  The  equation 
(2.8)  relating  the  frequency  u  and  the  wavenumber  n  of  propagating  waves  is  known 
as  the  dispersion  equation  for  the  medium.  In  the  case  of  vacuum  or  free  space  where 
9p  =  9c  =  0  so  that  n  =  co/c,  we  obtain  the  corresponding  phase  velocity  vp  =  c  =  the 
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speed  of  light  as  expected.  More  generally  the  phase  velocity  in  a  dielectric  medium 
with  conductivity  and  polarization  is  given  by 


vp  =  c/\Jl+igc/e0u  +  gp/e0.  (2.9) 

In  light  of  (2.9)  and  the  definition  of  a  dispersive  medium,  we  see  that  if  either 
gju  or  gp  depend  on  to,  we  will  have  dispersiveness.  Several  special  cases  are  worthy 
of  note. 

For  instantaneous  conductivity,  that  is,  gc(t,z)  =  a8(t)  so  that  (2.3)  reduces  to 
Ohm’s  Law  Jc  =  aE,  we  see  that  the  term  igc/e0co  becomes  io /eQu.  Thus  a  medium 
with  simple  Ohm’s  Law  conductivity  will  be  dispersive  (it  is  also  dissipative  in  the 
usual  sense  since  the  conductivity  term  in  (2.1)  becomes  ^E).  For  instantaneous 
polarization  (often  assumed  in  standard  treatments  of  the  Maxwell  theory)  we  find 
gp(t,  z )  =  eoX<5(t),  where  x  is  the  dielectric  susceptibility  constant  and  hence  gp/e o  =  x 
and  the  medium  is  not  dispersive.  One  must  turn  to  more  complicated  (and  more 
realistic)  models,  such  as  those  of  Debye  or  Lorentz,  to  have  a  polarization  based 
contribution  to  dispersiveness  in  a  medium.  For  the  usual  Debye  polarization  model 
[E1193,  p.386]  one  has 

gp(t)  =  e~t/Te0(es  -  )/r,  t  >  0,  (2.10) 

where  r  is  a  relaxation  parameter  and  es,  are  familiar  dielectric  constants.  In  this 
case  one  finds 

gp{uj)  =  c0(es-eoo)[i1+^J2_ 

For  the  Lorentz  model  [RMC93,  p.496]  we  have 

e  co2 

gp(t)  =  ~^2t  sin^ot,  t  >  0,  (2.11) 
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where  vq  =  \Ju'q  —  1/4 r2.  In  the  frequency  domain  this  yields 

9p(u)  =  e0  u)p[ 


4  T2 


=  e0upiy0 


(1  +  i2ru)2  +  4t2u'q- 
uj'2)  —  iu/r 


H 


(u>g  -  LU2)2  +  U2/t2 


and  again  we  have  a  polarization  based  dispersive  medium.  Higher  order  (the  Debye 
and  Lorentz  models  correspond  to  first  and  second  order,  respectively,  differential 
equation  models  for  the  polarization  P  -  see  [BBLOO]  and  the  references  therein) 
models,  as  well  as  combinations  of  such  models  also  lead  to  dispersion  in  a  medium. 

Thus,  in  summary  we  see  that  instantaneous  conductivity  but  not  instantaneous 
polarization  yields  dispersiveness  in  a  medium.  For  a  polarization  contribution  to  dis¬ 
persiveness  one  must  include  first  or  higher  order  polarization  models  (instantaneous 
polarization  can  be  correctly  viewed  as  zero  order  polarization  dynamics).  For  our 
semigroup  presentation  in  the  next  section  we  shall  therefore  consider  the  model  (2.1) 
with  instantaneous  conductivity  and  a  general  (higher  order)  polarization  model  given 
by  (2.2)  with  gp  =  g  where  the  DRF  g  is  assumed  smooth  in  time  (i.e. ,  without  loss 
of  generality  we  can  assume  no  instantaneous  component  for  g).  Such  distributed  pa¬ 
rameter  systems  are  of  interest  since  they  are  neither  simple  hyperbolic  nor  parabolic 
in  nature. 

For  simple  Ohm’s  Law  conductivity  and  instantaneous  polarization  (or  no  polar¬ 
ization),  the  system  (2.1)  becomes  a  well  understood  dissipative  or  damped  hyperbolic 
system  for  which  a  semigroup  formulation  can  readily  be  found  in  the  research  liter¬ 
ature  on  distributed  parameter  systems.  However,  for  (2.1)  with  polarization  based 
dispersiveness,  we  obtain  a  system  with  behavior  of  solutions  that  are  neither  stan¬ 
dard  hyperbolic  (finite  speed  of  wave  propagation  along  characteristics)  nor  standard 
parabolic  (infinite  speed  of  propagation  of  disturbances).  Indeed  for  (2.1)  with  either 
Debye  or  Lorentz  polarization,  rather  fascinating  solutions  can  be  observed.  These  in- 
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volve  the  formulation  of  so-called  Brillouin  and  Sommerfeld  precursors  where  a  pulsed 
excitation  (containing  waves  with  a  range  of  frequencies)  evolves  into  waves  propa¬ 
gated  with  different  velocities  which  coalesce  into  wave  “packets”  (see  Chapter  4  of 
[BBLOO]  and  [APM89]  and  the  references  therein  for  discussions  of  these  phenomena). 

It  is  of  both  mathematical  and  practical  interest  to  know  whether  these  interesting 
systems  can  be  described  in  a  semigroup  context.  The  potential  advantages  afforded 
by  a  semigroup  formulation  are  widespread  since  there  is  a  tremendous  literature  for 
control,  estimation  and  identification,  and  stabilization  of  systems  in  a  semigroup 
setting.  Results  for  both  stochastic  and  deterministic  control  methodologies  (in  both 
time  domain  and  frequency  domain)  including  open  loop  and  feedback  formulations 
are  abundant  [CZ95],  [BDDM92],  [BDDM93],  [CP78],  [vK93]. 

In  the  next  section  we  present  a  semigroup  formulation  of  the  system  (2.1)  with 
simple  Ohm’s  Law  conductivity  along  with  general  polarization  based  dispersiveness 
generated  by  polarization  laws  of  the  form  (2.2).  To  be  more  precise,  we  take  (2.1)  for 
t  >  0  and  z  G  (0, 1)  with  Jc(t,  z )  =  a(z)E(t,  z)  where  a(z)  vanishes  outside  C  (0, 1]. 
The  closed  region  fl  is  some  dielectric  material  region  (e.g.,  a  slab  or  several  slab-like 
regions)  containing  instantaneous  conductivity  as  well  as  non  trivial  polarization  of 
the  form  (2.2)  with  gp(t,z)  =  g(t,z )  vanishing  outside  z  G  Q.  Using  this  form  of 
conductivity  and  polarization  in  (2.1),  we  obtain  the  system 

E(t,  z)  +  —(a(z)  +  ^(O,  z))E(t,  z)  +  —0(0,  z)E(t,  z) 

Co  Co 

ft  l  l  . 

+  /  — g(t  —  s,  z)E(s,  z)ds  —  c2E"(t,  z)  = - Js(t,z).  (2.12) 

JO  e0  Co 

With  this  system  we  take  boundary  conditions  (see  [BBLOO]  for  details)  that  represent 
a  total  absorbing  boundary  at  z  =  0  and  a  supraconductive  boundary  at  z  =  1.  This 
can  be  expressed  by 

E(t,0)-cE'(t,0)  =  0  (2.13) 
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E(t,  1)  =  0. 


(2.14) 


With  the  definitions 

a(t,z)  =  —g(t,  z),  0(z)  =  —g(O,z) 

Co  Co 

7  (z)  =  —(°(z)  +  g(Q,z)),  J(t,z)  =  ~—js{t,z), 
e0  e0 

we  can  write  equation  (2.12)  as 

E  +  7 E  +  (3E  +  a  *  E  —  c2E"  =  J  (2.15) 

where  a*  E  is  the  usual  convolution 

a*E(t,  z)  =  [  a(t,  —  s,  z)E(s,  z)ds.  (2.16) 

Jo 

One  can  use  the  boundary  conditions  (2.13)  -  (2.14)  to  write  (2.15)  in  weak  or 
variational  form  so  as  to  seek  solutions  t  — >  E(t )  in  V  =  -##(0, 1)  =  {<fi  £  ^1(0, 1)  : 
4>(1)  =  0}  in  a  Gelfand  triple  setting  V  ^  H  V*  with  pivot  space  H  =  L'2( 0, 1). 
Under  modest  regularity  assumptions  on  a,  /3, 7  and  J ,  one  can  establish  existence, 
uniqueness  and  continuous  dependence  (on  initial  conditions  and  input)  of  solutions. 
Details  are  given  in  Chapter  3  of  [BBL00]. 
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3  A  Semigroup  Formulation 


We  turn  in  this  section  to  a  semigroup  formulation  for  the  dispersive  system  (2.12) 
-  (2.14)  or  equivalently,  (2.13)  -  (2.15),  with  instantaneous  conductivity  and  general 
(non  instantaneous)  polarization. 

For  our  development  we  assume  that  7, /3  G  L°°(0, 1)  while  a  G  L°°((0,T)  x 
(0, 1))  and  a,/3, 7  vanish  outside  fl  We  moreover  assume  that  a  can  be  written  as 
a(t,z )  =  ai(t)a2(z)  where  0  <  07  <  a2(z)  <  ctu  onOc  (0, 1]  for  positive  constants 
07,  au,  with  a2  vanishing  outside  fl  We  assume  that  t  — aq  (t)  is  positive,  monotone 
nonincreasing,  and  in  i/1(0,T)  so  that  aq (t)  <  0.  This  monotonicity  assumption 
is  typical  of  the  usual  assumptions  in  displacement  susceptibility  kernels  (e.g.,  see 
[BI08I,  p.  102]  or  [Hop77]).  We  shall  return  to  discuss  this  monotonicity  requirement 
further  after  our  semigroup  arguments  of  this  section. 

We  consider  the  term  (2.16)  given  by 

[  a(t  —  s)E(s)ds  =  [  a(t  —  s)E(s)ds 

JO  2—oo 

from  (2.15)  and  note  that  it  can  be  equivalently  written 

/  a(t  —  s)E(s)ds  =  f  a(—^)E(t  +  ^)d^ 

J  —  OO  J  —  OO 

«  f  <*(-z)E(t + m  =  f  G(0E(t + m 

J —r  J  — r 

where  G(£)  =  «(— £).  We  denote  Gi(^)  =  Qq(— ^)  so  that  G(£)  =  Gi(^)a2- 

The  approximation  is  valid  for  r  sufficiently  large  (r  =  00  is  permitted)  so  that 
a(t)  ~  0  for  t  >  v .  We  observe  at  this  point  that  Gi(£)  >  0  with  Gi(^)  >  0  on  (— r,  0]. 

As  introduced  in  the  previous  section,  we  take  V  =  72^(0, 1),  H  =  L'2( 0, 1)  with 
V  ^  H  ^  V*.  We  shall  have  use  of  H  =  L^2(Q),  the  space  L2(Q)  with  weighting- 
function  a2,  which  is  readily  seen  to  be  equivalent  to  E2(Q)  due  to  the  upper  and  lower 
bounds  on  a2  G  L°°(Q).  We  shall  denote  the  restriction  of  functions  (f>  in  L2(0, 1)  to 
Q  again  by  (j)  and  write  <fi  G  L2(fi)  or  <j>  G  L2  (fi)  whenever  no  confusion  will  result. 
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Using  the  above  definitions  and  approximating,  we  may  write  (2.15)  as 

E(t)  +  jE(t)  +  0E(t )  +  G(£)E(t  +  -  c2E"{t )  =  J(t).  (3.1) 

J—r 

Using  an  approach  given  in  [BFW88],  [BFW89],  [FI90]  and  [BMZ96]  for  viscoelastic 
systems,  we  define  an  auxiliary  variable  w(t)  in  W  =  L2Gi  (— r,  0;  H)  by  w(t)(9)  = 
E(t)  —  E(t  +  9),  —r<9<  0.  Since  G(9,  z)  >  0  for  9  G  (— r,  0],  z  G  52  we  may  take  as 
an  inner  product  for  W  the  weighted  L 2  inner  product 

(V,w)w=  [  Gi(9){r](9),w(9))jjd9  =  [  G?i(0)  f  a2{z)r](9,  z)w(9,  z)dz 

J  —r  J  —r  Jfl 

under  which  W  is  a  Hilbert  Space.  We  note  that  by  our  notational  convention  ex¬ 
plained  above,  we  have  w(t)  £  W  for  any  E(t,z )  with  E( -,  •)  e  LG  (—r,  0;  H ).  Using 
a  standard  shift  notation,  we  may  write  w(t)  =  E(t)  —  E(t  +  9)  =  E(t)  —  Et(9)  where 
E t{9)  =  E(t  +  9)  for  —  r  <  9  <  0.  Adding  and  subtracting  appropriate  terms  in  (3.1), 
we  find 

E(t)  +  ^E(t)  +  (3E(t)  +  f  G(0E(t)dti  -  f°  G(0[E(t)  -  E\0}^  ~  <?E"(t )  =  J(t) 

J —r  J —r 

or,  equivalently 

E(t)  +  7E  +  (0  +  Gu)E(t)  -  f°  G(0w(m)d£  -  c2E"(t)  =  J(t)  (3.2) 

where  Gn(z)  =  J°_rG(^)d^  =  a2(z)  /°rGi(0<^  and  w(t)(£)  =  E(t)  -  E t(£).  We 
observe  that  Gn,  like  p.  is  in  L2(Q)  and  L2( 0, 1). 

For  our  semigroup  formulation,  we  consider  (3.2)  in  the  state  space  Z  —  V  x  H  x 
IF  =  if^(0, 1)  x  L2(0, 1)  x  LGi(— r,  0,H)  with  states  (<; b,ip,r /)  =  (E(t),  E(t),w(t))  = 
(E (t) ,  E (t) ,  E (t)  —  E *(■)).  To  define  an  infinitesimal  generator,  we  begin  by  defining 
a  fundamental  set  of  component  operators.  Let  A  (E  G(U,  V*)  be  defined  by 

A(p  =  c2  <f>"  —  (/3  +  Gu)(j)  +  c2(j)'(0)5o  (3.3) 
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where  (50  is  the  Dirac  operator  =  ip( 0).  Then  we  find 

{  -  MA)v*,v  =  (CV,^%  +  {{P  +  Gu)4>,^)h  (3.4) 

so  that  it  is  readily  seen  that  d\  :  V  x  V  t— )•  C  defined  by 

^i(0,VO  =  (  -  M^)v\v  (3-5) 

is  symmetric,  V  continuous  and  V  coercive  (i.e. ,  <7i(^>,</>)  >  Ci\<j)\y  —  M&Ih  for  con¬ 
stants  Ao  and  c\  >0). 

We  also  define  operators  B  G  C{V,  V *)  and  K  G  £(W,  if)  by 

B(fi  =  —7^  —  c<^(0)^o  (3.6) 


so  that 

(  -  +  c<p(0)ip(0) 

and,  for  7]  G  W  =  EGi(—r ,  0;  if), 

0  z  g  [0,  i]\n 

,  S-r  G(f) v(0<%  *  6  a 
Since  &’(£.  z)  =  0  for  z  G  [0, l]\fl,  we  abuse  notation  and  write  this  as 


(Kt/)(z)  =  { 


Kr,  =  f  G(Z)ri(m 
J—r 


(3.7) 


(3,8) 


even  though,  strictly  speaking,  r)(£,  z )  is  only  defined  for  zGfi. 

With  these  definitions  and  notations,  equation  (3.2)  can  then  be  written  as 

(E,  <f>)y,  y  +  (  —  AE,  4>)y*  y  +  (  —  BE ,  < t>)y*  y 

+  (-K(E-Et),4>)yty=(j,(f>)yty 


or 

E{t)  =  AE(t)  +  BE(t)  +  K(E(t)  —  Ef)  +  J(t)  in  V*.  (3.9) 
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We  rewrite  equation  (3.9)  as  a  first  order  system  in  the  state  £(f)  =  (. E(t ),  E(t),w(t )) 
where  w(t)  =  E(t)  —  El .  To  aid  in  this  we  introduce  another  operator 

D  :  dom  D  C  W  W 

defined  on  dom  D  =  {q  G  if1(— r,  0;  H)\rj(0)  =  0}  by 

DV{0)  =  j^(0)- 

We  then  observe  that  u>(t)  =  E(t)  —  Et  satisfies 

w(t)(0)  =  E(t )  -  £(t  +  0)  =  £(t)  -  DE*(0) 

=  E(t)  +  D(£(f)  -  Et(9))  =  E(t)  +  Dw(t)(0). 

Thus  we  may  formally  rewrite  (3.9)  as  a  first  order  system  and  adjoin  to  it  the 
equation 

w(t)  =  Dw(t )  +  E(t).  (3.10) 

We  then  obtain  the  first  order  system  for  £(t)  given  by 

C  (t)=AC(t)  +  E(t)  (3.11) 


where  A  given  by 


is  defined  on 


0/0 
A=  A  B  K 
L  0  ID 


(3.12) 


dom  A  —  {((/),  ip,  rj)  G  Z\ip  G  V,rj  G  dom  D ,  Acp  +  Bip  G  H}.  (3.13) 


That  is,  A&  =  (<p,  A(p  +  Bip  +  Kr /,  ip  +  Drj)  for  $  =  (<p,  ip,  rf)  in  dom  A.  The  forcing 
function  E  in  (3.11)  is  given  by  E  =  col(0,  J ,  0).  To  argue  that  A  is  the  infinitesimal 
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generator  of  a  Co-semigroup,  we  actually  consider  the  system  (3.11)  on  an  equivlaent 
space  Z\  =  V\  x  H  x  W  where  V\  is  the  space  V  with  equivalent  inner  product 

,  <^2) vy  =  &i(<Pi,<fa)  where  d\  is  the  sesquilinear  form  given  in  (3.5).  Recall  that 
<7i  is  symmetric,  V  continuous  and  V  coercive  so  that  it  is  topologically  equivalent 
to  the  V  inner  product. 

We  are  now  ready  to  prove  the  following  generation  theorem. 

Theorem  Suppose  that  7,  [i  £  L°°(0, 1),  a  (E  L°°((0, 1)  x  (0, 1))  with  a,  (1,  7  vanishing 
outside  fh  We  further  assume  that  a  can  be  written  a(t,z )  =  a\(t)a2(z)  where 
aq  (E  H\0,T)  with  aq(t)  >  0,  aq(t)  <  0,  and  0  <  07  <  a2(z)  <  au  for  positive 
constants  aL,au-  Then  the  operator  A  defined  by  (3.12),  (3.13)  is  the  infinitesimal 
generator  of  a  Co-semigroup  on  Z\  and  hence  on  the  equivalent  space  Z. 

To  prove  this  theorem,  we  use  the  Lumer-Phillips  theorem  ([Paz83,  p.14]) .  Since 
Z\  is  a  Hilbert  space,  it  suffices  to  argue  that  for  some  Ao,  A  —  Ao /  is  dissipative  in 
Zi  and  1Z(XI  —  A)  =  Z\  for  some  A  >  0,  where  V,(X  1  —  A)  is  the  range  of  XI  —  A. 
We  first  argue  dissipativeness. 

Let  $  =  ((f),  ip,  rj)  (E  dom  A.  Then 

(A®^)Zl  =  <^M)vi  +  (A(P  +  B*p  +  Kg,ip)H  +  (xp  +  Dg,g)w 

=  (VN  <P)Vl  +  (M  +  BxPi  $)v,v  +  (Kth  *P)h  +  +  DSi  v)w 

=  °i (ip,  <t>)  ~  &i(<P, i>)  +  (Bip,  ip)H  +  (Kg,  ip)H  +  (ip  +  Dg,  g)w 
=  -('Y'fp,  tp)H-  cl^(°)  I2  +  (KV,  'tP)H  +  ('lP  +  DV,  V)w 
<  ItIooIV’IIt \(Kg,ip)H\  +  \{ip  +  Dg,g)w\.  (3.14) 

We  consider  estimates  for  the  last  two  terms  in  (3.14)  separately.  From  (3.8)  we 
have 

\(Kg,  ip)H\  =  |  f  G1(d)(g(d),'tp)gd9\ 

J  —r 
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<  r  omnm^e 

J  —r 

<  \f_G,mW)W  +  W\)de 

<  h\v\w  +  k^\2H- 

Moreover, 

\(*l>,v)w\  -  [  Gi (0)\(fav(0))H\d0 

J  —r 

<  /_°  a.mlm+lmmdf) 

<  fa\i>\2H  +  h\v\w- 

Finally,  since  G\{9)  >  0,  G\  >  0,  and  rj  G  dom  D  requires  77(0)  =  0,  we  may  argue 

(DV,V)W  =  f  G1{6){Dr]{e),r]{e))ilde 

J  —r 

=  Ca'Kl  Me^de 

=  f_r^Gm\m\\)de  -  f_r\(GMm\\)M 

=  G1(0)1|,(0)||  -  G,(— r)i|,(— r)| \  -  1  faMmilde 

<  0. 

Combining  these  estimates  with  (3.14),  we  obtain  for  $  G  dom  A 

(A®,  $)Zi  <  hlocM  +  +  k4)\ri\2w  +  (k2  +  k3)\Tp\l 

<  Aoj<h||i, 

which  yields  the  desired  dissipativeness  in  Z\. 

To  establish  the  range  statement,  we  must  argue  there  exists  some  A  >  0  such 
that  for  any  given  T  =  (/q  zv,  £)  in  Z,  there  exists  <f>  in  dom  A  satisfying 

(A/-^)$  =  T.  (3.15) 
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In  view  of  the  definition  of  A,  the  equation  (3.15)  is  equivalent  to  the  system 

A  (f>  —  ip  =  n 

-A(j>  +  (A  —  B)xp  —  Kt]  =  v  (3.16) 

-V>+(A -D)r]  =  t 

for  (cf),  if),  rf)  G  dom  A,  (/i,  v,  £)  G  Z  =  VxHxW.  The  first  equation  is  the  same  as  ip  = 
X(f)  —  n  while  the  third  can  be  written  as  77  =  (A  — D)_1(^  +  ^)  =  (A  —  D)-1(f  +  X(p  —  //). 
These  two  equations  can  be  substituted  in  the  second  to  obtain  an  equation  for  cf>.  If 
this  equation  can  be  solved  for  <f>  G  V,  then  the  first  and  third  can  then  be  solved  for 
if)  and  r],  respectively.  The  equation  for  <f>  that  must  be  solved  is  given  by 

-A<t>  +  (A  -  B)(\<j>  -id)-  K(  A  -  D)~l(i  +  A  <f>-n)  =  v 

or 

[A2  -  A  B-A-  K(  A  -  D)_1A]  (f)  =  (A  -  B)/i  +  v  +  K(  A  -  D)-1(£  -  /i).  (3.17) 

If  we  can  invert  (3.17)  for  cf)  G  V,  then  ip  =  A0  —  /d  is  in  V,  77  =  (A  —  D)-1  [£  +  A</i>  —  /x] 
is  in  dom  D  C  IT  and 

I^  +  S-0  =  A2</i-  AJn-^-/?(A-D)-1(^  +  A</)-/i) 

=  Xip  —  d  —  krj 

is  in  H  so  that  ((f),  ip,  T])  is  in  dom  A  and  solves  (3.16). 

Thus  the  range  statement  reduces  to  solving  (3.17)  for  (p  G  V .  This  in  turn  reduces 
to  invertability  of  the  operator  A2  —  A B  —  A  —  k(X  —  D)~x A. 

We  first  observe  that  (A  —  D)-1  =  (1  —  ex')/X  since  (A  —  D)(  1  —  exd )  =  X  while 
r](0)  =  [£  +  A (p  —  n]  satisfies  77(0)  =  0  and  hence  is  in  dom  D. 

Thus,  for  <p  G  H,  K( X  —  D)~l X  satisfies 

(k(l-exe)<P,<P)H  =  f  Gl(e)(l-exe)(<P,<P)flde 

J  —r 

<  k5\(p\2H 
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and 


((A2  -  XB)(j),  4>)v»y  =  ((A2  +  X'l)^  +  Ac^(0)50,  4>)y»y 

=  <(A2  +  A7)0,  <f>)H  +  Ac| 0(0) |2 
>  h\<f>?H  for  A  sufficiently  large. 

Hence  for  A  sufficiently  large  we  have 

((A2  -  A  B- A-  I<(  A  -  D)~1X)(j),  </>)y,y 

=  ((A2  -  A 4>)v,y  +  &M,  4>)  -  (K( A  -  D)~1\(j),  <j>)H 

>  ke\(f>\2H  +  Ci\(f>\y  —  Xo\(f>\2H  —  ks\(f)\2H 

=  Ci \(j)\y  +  (A;6  —  Ao  —  ks)\(j)\2H. 

Thus  if  we  define  the  sesquilinear  form 

Oxtt, i>)  =  ((A2  —  XB  —  A  —  K( A  -  D)~1\)4>,  ^)v.jV, 

we  see  that  for  A  sufficiently  large,  cr*  is  V  coercive  and  hence,  by  the  Lax-Milgram 
lemma  [Wlo87],  it  is  invertible.  It  follows  immediately  that  (3.17)  is  invertible  for 
4>  6  V.  This  completes  the  arguments  to  prove  the  Theorem. 

Let  S(t)  denote  the  semigroup  generated  by  A  so  that  solutions  to  (3.11)  are  given 
by 

C(t)  =  S(t)Co+  [  S(t  —  s)!F(s)ds.  (3.18) 

Jo 

Solutions  are  clearly  continuously  dependent  on  initial  data  Co  and  the  nonhomoge- 
neous  perturbation  T.  The  first  component  of  £(t)  is  a  generalized  solution  E(t) 
of  (3.1).  One  can  now  argue  that  the  solution  agrees  with  the  unique  weak  so¬ 
lution  obtained  in  Chapter  3  of  [BBL00],  by  using  the  arguments  in  Chapter  4.4 
of  [BSW96].  Briefly,  one  argues  equivalence  for  sufficiently  regular  initial  data  and 
nonhomogeneous  perturbation.  Then  density  along  with  continuous  dependence  is 
used  to  extend  the  equivalence  to  more  general  data  (see  [BSW96]  for  details). 
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4  Concluding  Remarks 


In  the  previous  section  we  presented  a  semigroup  generation  theorem  under  general 
conditions  on  the  coefficients  at,  (3, 7  of  (2.15).  The  only  possibly  restrictive  condition 
involved  a(t.  z )  =  j^g(t,  z )  =  a1(t)a2(z)  where  it  is  required  that  aq(t)  >  0,  di(t)  <  0. 
We  consider  more  closely  the  condition  for  some  common  polarization  laws. 

For  Debye  polarization  in  a  region  fi,  we  have  cti(f)  =  j^gP(t)  where  gp  is  given  in 
(2.10).  That  is, 

gp(t)  =  e0 - - - e  > 

so  that 

ati (t)  =  —  gp(t)  =  —  e°°K~t/r  >  0 
e0  ts 

and 

dl(t)  =  .fo-  C°°)e-«/T  <  0. 

T4 

Thus  Debye  polarization  satisfies  the  conditions  of  the  generation  theorem  and  the 
associated  system  generates  a  Co  semigroup. 

For  Lorentz  polarization,  we  have  (recall  (2.11)) 

gp(t)  =  ^^e“t/2rsin  vot 
vo 

and  hence 

«i  (t)  =  —  gp(t)  =  —  [(77W  -  Wq)  sin  ^0 1  ~  —  cos  vot] . 

We  therefore  see  that  it  is  not  possible  to  conclude  that  a\{t )  >  0  or  di(t)  <  0  so  that 
our  generation  theorem  does  not  guarantee  a  semigroup  representation  for  systems 
with  a  Lorentz  polarization  law.  In  spite  of  this,  we  do  believe  that  the  Lorentz 
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law  does  yield  a  system  with  a  semigroup  representation.  We  conjecture  that  the 
proof  of  the  theorem  we  present  can  be  modified  to  weaken  the  hypothesis  on  a  so 
as  to  include  Lorentz  and  other  oscillatory  (even  order)  polarization  models.  We  are 
currently  pursuing  these  ideas. 

In  closing  we  point  out  that  the  general  class  of  dielectric  response  functions 
consisting  of  a  linear  combination  of  decreasing  exponentials  (essentially  multiple 
Debye  mechanisms)  suggested  for  glasseous  materials  by  Hopkinson  [Hop77]  (see  the 
discussion  in  [BI08I,  p.  101-103])  are  included  under  the  theory  presented  in  this  note. 
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